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Abstract. Since the launch of the Large Area Telescope (LAT) on board the Fermi spacecraft 
in June 2008, the number of observed 7-ray pulsars has increased dramatically. A large number 
of these are also observed at radio frequencies. Constraints on the viewing geometries of 5 
of 6 7-ray pulsars exhibiting single-peaked 7-ray profiles were derived using high-quality radio 
polarization data [1]. We obtain independent constraints on the viewing geometries of 6 by 
using a geometric emission code to model the Fermi LAT and radio light curves (LCs). We find 
fits for the magnetic inclination and observer angles by searching the solution space by eye. Our 
results are generally consistent with those previously obtained [1], although we do find small 
differences in some cases. We will indicate how the 7-ray and radio pulse shapes as well as their 
relative phase lags lead to constraints in the solution space. Values for the flux correction factor 
(/n) corresponding to the fits are also derived (with errors). 


1. Introduction 

We study 6 pulsars (J0631+1036, J0659+1414, J0742-2822, J1420-6048, J1509-5850, and 
J1718— 3825) detected by Fermi LAT, all exhibiting single-peaked pulse profiles within current 
statistics at energies >0.1 GeV. These pulsars are also detected in the radio band. Both of these 
properties aid in constraining the possible solution space for the respective pulsar geometries, 
as derived from predictions of their light curves (LCs). 

A study of the pulsars’ geometric parameters, the inclination and observer angles a and 7 
has been performed [1] using the radio polarization and LC data. They constrained the solution 
space for a and /3 of 5 of the individual pulsars, with the impact angle (3 — ( — a, using fits to 
these radio data as well as predictions for the value of the half opening angle, p , of the radio 
beam derived from the radio pulse width (e.g., [ 2 ]). 

The aim of this study is to obtain similar constraints on a and £ of the same pulsars using 
an independent, multiwavelength approach. We use a geometric pulsar emission code (e.g., [3]) 
to model the Fermi LAT and radio LCs, and fit the predicted radio and 7 -ray profiles to the 
data concurrently, thereby significantly constraining a and £. This also allows us to test the 
consistency of the various approaches used to infer the pulsar geometry. 


2. Model 

We use an idealized picture of the pulsar system, wherein the magnetic field has a retarded 
dipole structure [4] and the 7 -ray emission originates in regions of the magnetosphere (referred 
to as ‘gaps’) where the local charge density is sufficiently lower than the Goldreich-Julian charge 
density [5]. These gaps facilitate particle acceleration and radiation. We assume that there are 
constant-emissivity emission layers embedded within the gaps in the pulsar’s corotating frame. 
The location and geometry of these emission layers determine the shape of the 7 -ray LCs, and 
there exist multiple models for the geometry of the magnetosphere describing different possible 
gap configurations. 

We included two models for the 7 -ray emission regions in this study, namely the Outer 
Gap (OG) and Two-Pole Caustic (TPC) models. In both the OG and TPC models (the Slot 
Gap model [ 6 ] may serve as physical basis for the latter), emission is produced by accelerated 
charged particles moving within narrow gaps along the last open magnetic field lines. In the 
OG model [7], radiation originates above the null charge surface (where the Goldreich-Julian 
charge density changes sign) and interior to the last open magnetic field lines. The TPC gap 
starts at the stellar surface and extends along the last open field lines up to near the light 
cylinder 1 , where the corotation speed approaches the speed of light [ 8 ] . The special relativistic 
effects of aberration and time-of-ffight delays, which become important in regions far from the 
stellar surface (especially near the light cylinder), together with the curvature of the magnetic 
field lines, cause the radiation to form caustics (accumulated emission in narrow phase bands). 
These caustics are detected as peaks in the observed 7 -ray LC [ 8 , 9]. 

We used an empirical radio cone model [10], where the beam diameter, width, and altitude 
are functions of the pulsar period P, its time derivative P, and the frequency of observation v 1 
in order to obtain predictions for the radio LCs. This is different from the approach taken by 
[ 1 ], as we have a different prescription for the radio emission altitude and p, and we do not use 
the polarization data. 

3. Method 

Using the above geometric models, we generated LCs as a function of a and £, keeping the gap 
widths and extents constant. Due to the large size of the (a,£)-space at a 1° resolution, it is 
impractical to search blindly for the optimum LC by eye. We therefore started by generating a 
so-called atlas of LCs with a 10° resolution for each pulsar and first identified candidate LC fits. 
We then performed refined searches in the regions of (a,£)-space around those candidate fits at 
a 5° and later 1° resolution. This approach of generating all the LCs in a candidate region and 
comparing them directly, allowed us to obtain constraints on a and We then inferred contours 
for a and () for each pulsar and for each geometric model. We lastly used these contours to find 
the corresponding values of fn (with errors; see Section 3.2) predicted by the OG and TPC 
models. 

3.1. Obtaining the contours 

To obtain an (a,£)-contour delineating the region of plausible LC fits, we start at the candidate 
solution (5° resolution) and first move away from it in steps of 1° in a. At each step we then 
step through £ until the upper and lower bounds of the contour are located for the fixed a. This 
procedure is repeated until no further acceptable LCs are found. As an example, the red regions 
A and B in the left panel of Figure 1 indicate the resulting contours obtained in this manner 
for PSR J1509— 5850 for the TPC case. 

1 The slight difference in transverse polar position of the gaps with respect to the magnetic axis (as motivated by 
physical models) will result in a systematic shift in best-fit (a,£)-contours between the two models. See Section 4. 
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Figure 1 . Left: The (a,£)-contours (regions A and B) for radio and 7 -ray LC fits of PSR 
J1509— 5850 for the TPC model (red). The dark grey and light grey regions correspond to 
double-peaked and single-peaked radio LCs respectively. The solid lines indicate approximately 
the locations of the boundaries to the contours due to: (a) the relative phase lag between the 
radio and 7 -ray peaks being too large, ( 6 ) the bridge emission’s relative intensity being too 
low, and (c,d) the 7 -ray peak being too narrow. The dashed line (bounded area) indicates the 
resulting errors on a and with the green block indicating our best fit. Right: Representative 
LC fits corresponding to the contours A and B. 


Oftentimes we find two disjoint solution contours, with one of the two sometimes producing 
a better LC fit. The example shown here is one of these cases, as can be seen from the 
representative LCs shown in Figure 1. Contour B yields a considerably better fit than contour A, 
and thus we can ignore contour A in favour of contour B as indicated by the dotted line in 
Figure 1. The values of a and f reported in Table 1 are those at the centre of this dotted box, 
while the errors on these values are conservatively chosen to encompass the dimensions of the 
box. In this case the resulting values are a — (61° ± 5°) and f — (44° ± 7°). 


3.2. Finding fa 

It is important to be able to convert the observed energy flux of a pulsar to its all-sky luminosity. 
The flux correction factor Jq is used for this purpose. It is a highly model-dependent parameter, 
and allows us to determine what fraction of the all-sky luminosity we are observing from a pulsar 
if the geometry is known (a and £), for a given model. It is also crucial in deriving the 7 -ray 
efficiency of a pulsar, and is defined as [ 11 ] 


/n(«, Ce) 


Jf F 7 (a, (,f>) sin (d(d<f> 
2fF 7 (a,CE,<f>) dcf> 


(1) 


with F 7 being the photon flux per solid angle (‘intensity’), and Qe the Earth line-of-sight. The 
value of /h is typically taken to be 1 , meaning that the observed energy flux is assumed to be 
equal to the average energy flux over the entire sky. Due to the OG model predicting emission 
over a relatively smaller region of phase space, its associated fn values are typically smaller 
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Figure 2. Contour plot of /q as function of a and £ implied by the TPC model for PSR 
J1509-5850. The red area indicates the best-fit (a, £)-contour for the radio and 7 -ray LCs. 


than those of the TPC model for the same a and (Note, we assume that F 7 (£, <j>)/F^ t tot is 
approximately equal to the ratio of observed energy flux vs. total energy flux.) 

We can now use the (a , ^-contours obtained by the method discussed in Section 3.1 to 
constrain the value of /n by computing fa(a, Q and overplotting those contours (see Figure 2). 

4. Results 

We infer values for a and ( for each pulsar with typical errors of ~ 5° (see Table 1). The 
tabulated values are the average values of a, £, and j3 implied by the solution contours, while 
the errors are chosen conservatively so as to include the full (non-rectangular) contour. We note 
that once a solution is found at a particular a and £), it is worthwhile to study the LC at the 
position a! — C, and (,' — a (i.e. , when the angles are interchanged). Such a complementary 
solution may provide a good fit in some cases (see the alternative solutions listed in Table 1), 
due to the symmetry of the model as well the symmetry in pulsar geometry (i.e., |/3| remaining 
constant under this transformation). 

As mentioned in Section 2, we see a shift in best-fit contours, with the TPC model consistently 
yielding smaller a and £ values compared to those obtained using the OG model due to a different 
gap geometry. This is because similar phaseplots are obtained when a larger a is chosen for the 
OG model than for the TPC model. Reproduction of the radio pulse shape then necessitates a 
larger value of £ for the OG fit, since the contour is constrained to appear in one of the grey 
regions (as determined by peak multiplicity; see Figure 1). The net effect of this is that the 
OG contours are found farther from the origin at a — £ — 0 along the grey bands than the 
corresponding TPC contours. 

Our results for /q for each of the pulsars (and models) are given in Table 1. 



5. Discussion and Conclusions 

The good constraints on a and £, obtained using the method described above, emphasizes the 
merit of a multiwavelength approach. For example, the light grey regions in Figure 1 indicate 
single-peaked radio profile LC solutions. These would represent typical constraints one would be 
able to derive for a and £ when considering only the radio profile shapes. Similarly, when only 
considering the (single-peaked) 7 -ray LCs, one would find relatively large (a, ^-contours [11]. 
The requirement of fitting both the radio and 7 -ray profile shape, as well as their relative phase 
lag, results in much smaller contours, as shown in Figure 1. 

Generally, our best-fit (a,£) compare favourably with those inferred by [1] for the first three 
pulsars (see Table 1). However, comparison is hampered by uncertainties in estimating the half 
opening angle, p, which sensitively influences the optimal solutions obtained by [1] . Even a small 
error of 5° on p leads to relatively large errors on the allowed a, so that our best fits would then 
be included in their inferred parameter ranges. Direct comparison of our results with those of 
[1] for PSR J1420— 6048 and PSR J1718— 3825 is not feasible due to different approaches for 
estimating the radio pulse width, W. We only model the most prominent radio peak, while [1] 
measures W as the width of the total two-peaked profile. Implementing a factor 2 difference in W 
will change the p-contours and thus their best-fit solution considerably, improving the agreement 
between the two approaches. Comparison with [1] is not possible for PSR J1509-5850 as the 
lack of polarization data inhibited inference of a best-fit (a,/3) by [1] in this case. 

The next step for this kind of study is to apply a mathematically rigorous method of 
determining the best-fit LC solutions as a function of many free model parameters to the six 
pulsars modelled in this paper. An example is the Markov chain Monte Carlo method which 
has been successfully applied to millisecond pulsar LCs [12]. 


Table 1. Values for a, £, /3, and /q derived for 6 pulsars for both the OG and TPC models. 
The *’s indicate alternate best-fit solutions upon interchange of a and 


Pulsar 

Model 

a (°) 

C(°) 

n°) 

fn 

J0631+1036 

OG 

74 ±5 

67 ±4 

- 6 ± 2 

0.93 ±0.06 


TPC 

71 ±6 

66 ±7 

- 5 ± 3 

1.04 ±0.04 

J0659+1414 

OG 

59 ±3 

48 ±3 

-12 ±5 

1.16 ±0.53 


TPC 

50 ±4 

39 ±4 

-13 ±6 

1.64 ±0.04 


TPC* 

38 ± 1 

50 ±4 

11 ± 4 

1.63 ±0.05 

J0742— 2822 

OG 

86 ±3 

71 ±5 

-16 ±6 

0.99 ±0.10 


OG* 

71 ±6 

86 ±4 

16 ±6 

0.81 ±0.09 


TPC 

64 ±8 

80 ±4 

15 ±6 

0.88 ±0.41 

J1420— 6048 

OG 

67 ±5 

45 ±7 

-22 ±9 

0.77 ±0.13 


TPC 

64 ±6 

43 ±8 

-21 ±9 

0.90 ±0.10 


TPC* 

42 ±5 

63 ±5 

21 ±9 

0.77 ±0.06 

J1509— 5850 

OG 

66 ±4 

50 ±7 

-18 ±8 

0.77 ±0.11 


TPC 

61 ±5 

44 ± 7 

-18 ±8 

0.89 ±0.10 

J1718— 3825 

OG 

67 ±6 

48 ±6 

-19 ±8 

0.76 ±0.12 


TPC 

61 ±5 

43 ±6 

-19 ±8 

0.86 ±0.07 


TPC* 

42 ±6 

62 ±5 

19 ±7 

0.83 ±0.10 
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